Quantum Phase Transitions in Capacitively Coupled Two-Dimensional 

Josephson-Junction Arrays 



00 
Oh- 

c/3 



Mahn-Soo Choi* 

Department of Physics and Pohang Superconductivity Center, Pohang University of Science and Technology, Pohang 790-784, 

Korea 
(September 4, 1998) 

Quantum-phase transitions in two layers of ultrasmall Josephson junctions, coupled capacitively 
with each other, are investigated. As the interlayer capacitance is increased, the system at zero 
temperature is found to exhibit an insulator-to-superconductor transition. It is shown that, unlike 
one-dimensional arrays with a similar coupling configuration, the transition cannot be accounted for 
exclusively by particle-hole pairs. 

PACS Numbers: 74.50.+r, 67.40.Db, 73.23.Hk 
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Capacitively coupled systems of charges have attracted 
significant attention in recent years, raising the possibil- 
ity of current drag effects: The current fed through either 
of the systems, owing to Coulomb interaction, induces a 
secondary current in the other system. Such a drag effect 
depends strongly on the dimensionality and the structure 
of the system in its mechanism and behavior. The cur- 
rent drag in two capacitively coupled two-dimensional 
(2D) electron gases |lj was attributed to momentum- 
transfer mechanism due to Coulomb scattering M and 
fairly small in magnitude. By contrast, recent theoretical 
prediction |3| and experimental demonstrations with 
two capacitively coupled one-dimensional (ID) arrays of 
submicron metallic tunnel junctions have shown that the 
primary and the secondary currents are comparable in 
magnitude but opposite in direction in a certain region 
of applied voltage. In such tunnel junction systems, the 
current drag is attributed to the transport of electron- 
hole pairs, which are bound by the electrostatic energy 
of the coupling capacitance. Lately, it has been suggested 
that the momentum-transfer mechanism can also lead to 
the absolute current drag in ID electron channels cou- 
pled electrostatically with each other ||. The current 
drag effects in capacitively coupled 2D arrays of tunnel 
junctions has not been studied and will be examined in 
this work. 

More interestingly, when the tunneling junctions are 
composed of ultrasmall superconducting grains, the 
counter part of the electron-hole pair becomes the pair 
of excess and deficit Cooper pair, which will be simply 
called the particle-hole pair. Furthermore, in such ul- 
trasmall Josephson junction systems the competition be- 
tween the charging energy and the Josephson coupling 
energy is well known to bring the novel effects of quan- 
tum fluctuations — 1]~0|] - Combined with these quantum 
fluctuation effects, the pair transport phenomena in cou- 
pled ID Josephson junction arrays (JJAs) has recently 
been proposed to drive an insulator-to-superconductor 
transition pi] ]. 

In this paper, two 2D arrays of ultrasmall Joseph- 



son junctions, coupled capacitively with each other, 
are considered. Quantum-phase transitions are exam- 
ined at zero temperature, focusing on the roles of the 
particle- hole pairs. The system is transformed into two 
three-dimensional (3D) systems of classical vortex loops, 
which are topologically coupled but otherwise indepen- 
dent of each other. The resulting model reveals that 
as the coupling capacitance increases, in appropriate re- 
gions of parameters, the system exhibits an insulator- 
to-superconductor transition. Contrary to the ID coun- 
terpart with a similar coupling scheme JlT| ], the transi- 
tion cannot be ascribed exclusively to the condensation of 
the particle-hole pairs. Accordingly, it is also remarked 
briefly that the accompanying drag of supercurrents in 
the superconducting phase is not absolute in general. In 
the vicinity of the transition, however, the particle-hole 
pairs still play major roles, and therefore current drag 
can be large. 

As a matter of fact, capacitively coupled 2D JJAs have 
been studied by several authors, but in different context 
and different regions of parameter space Besides, 
the capacitive coupling should be distinguished from the 
Josephson coupling such as considered in multi-layered 
systems The capacitively coupled JJAs can pre- 

sumably be realized in experiment by current techniques, 
which have already made it possible to fabricate sub- 
micron metallic junction arrays with large inter-array 
capacitances as well as large arrays of ultrasmall 
Josephson junctions 0. 

Each of the two arrays (£ — 1,2) of Josephson junc- 
tions considered here is characterized by the Josephson 
coupling energy Ej and the charging energies Eo=e 2 /2Cq 
and Ei=e 2 /2Ci, associated with the self-capacitance Co 
and the junction capacitance C\ , respectively (see Fig. [l]) . 
The two arrays are coupled with each other by the capaci- 
tance Cj, with which the electrostatic energy Ej=e 2 /2Cj 
is associated, while there is no Cooper-pair tunneling be- 
tween the arrays. The intra-array capacitances are as- 
sumed to be so small (Eq,Ei ^> Ej) that, without the 
coupling, both arrays would each be separately in the in- 



sulating phase m. It is also assumed that the coupling 
capacitance is sufficiently large compared with the intra- 
array capacitances; Cj Co, C\. In that case, the elec- 
trostatic energy of the particle-hole pair (~ Ei) is much 
smaller than that of an unpaired charge (~ E n ,Ei); the 
particle- hole pair, bound by the binding energy of order 
of Eq — Ei or E\ — Ej, is thus much favorable than the 
unpaired charges. For the most part, this work is devoted 
to the case of identical arrays, but non-identical arrays 
will also be briefly discussed. 

The system can be well described by the Hamiltonian 



H = -h E E n(£; r)c ~ 1{i > £ ' ; r ' r ' )n{£/; r ' } 

£,£' r,r' 

- 2K J2 E cos [</»(£; r)-^;r')], (1) 



where r = (x, y) denotes the 2D lattice vector in units of 
lattice constant, the coupling constant has been defined 
by 2K=y / E J /4E I , and the energy has been rescaled in 
units of the Josephson plasma frequency ftto p = ^AEjEj. 
The number n(£; r) of excess Cooper pairs and the phase 
4>(£; r) of the superconducting order parameter on the 
grain at r in array £ are quantum mechanically conju- 
gate variables: [n(£; r), <\>(l'\ r')]=iS Tr >5w ■ The capaci- 
tance matrix in Eq. (|IJ) takes the form 



(2) 



where the submatrices C(r, r') are defined by the Fourier 
transform C(q) = C +CiA(q) with A(q) = A(q x ) + 
A(gj,); A(fc) = 2(1— cos k). Here all the capacitances 
have been rescaled by the relevant capacitance scale 2Cj: 
C Q /2C I ^C and C 1 /2C I ^C 1 . 
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It is convenient to write the partition function of the 
system in the imaginary-time path integral representa- 
tion 

2=0 E / <Wr,r) exp[-S] (3) 

with the Euclidean action 

S= TrJ2 E n(^r,r)C- 1 (^^;r,r')n(^r',r) 



-2JfX;EE«»V J -^;r J r) 

l r,T j 

+i EE n ^ r ' r ) v -^' r ' r )' 



(4) 



where Vj (j = x, y) and V r denote the difference op- 
erators in the spatial and the imaginary-time directions, 
respectively, and the (imaginary-) time slice St has been 
chosen to be unity (in units of lo^ 1 ) [^5|. The highly 
symmetric form of Eq. (|]) with respect to space and time 
makes it useful to introduce the space-time 3-vector no- 
tation r = (r, r) , and analogous notations for all other 
vector variables. We then apply the Villain approxi- 



mation 1 16 1 to rewrite the cosine term as summation 
over integer fields {m x {£\ f), m y (I; ?)} = {m(£;r)}. Fur- 
ther, with the aid of Poisson resummation formula fig 1 ] 
and Gaussian integration, we rewrite the charging energy 
term as summation over another integer field {m T (£; r)}, 
to obtain the partition function 



Z . 



IIE 



;r) exp{-5} 



(5) 



with 



S = K C (^' e '> r > r ') S rr' i^r<f>(£; ?) - 2nm T (£; r)] [W T c/)(£'; r') - 2nm T (£'; ?)] 

£,£';r,r' 

+ K 4>(t; r)-27rm(£;r)] 2 . (6) 

£;r 

The variables </>(£; r) and m(£; r) can be usefully replaced by ^(r) = 0(1; r) ± </>(2; r) and m ± (f) = m(l; r) ±m(2; r), 
respectively. In this way, one decomposes the Euclidean action in Eq. (^|) into the sum S = S++S- with S ± defined 
by 

S ± = +^KY,C ± (r,r')S(T,T') [V T ^(f) - 27rm±(?)] [V T 4r t (f)-2irmf(f)] 

+ i/^[v^(f)-2,m ± (f)] 2 . (7) 
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Here the new capacitance matrices C ± (r,r') are defined 
via the Fourier transforms C + (q) = C(q) and C~(q) = 
1 + C(q), respectively. Now one follows the standard 
procedures to integrate out {0 ± (r)}. Then, apart 

from the irrelevant spin wave part, one can finally ob- 
tain the 3D system of classical vortex lines, which is also 
decomposed into two subsystems Hy — Hy 



H v with 



2 ^ K E E v t ( ? ) U±(?-?)v±(r>), (8) 



where the interactions between vortex line segments are 
defined via their Fourier transforms 



fif(<D = £^ 

11 A(q)+C±(q)A( W ) 

^(q) = 



A(q)+C±(q)AH 



(9) 
(10) 



Here the vortex lines are manifestation of the particle- 
hole pairs, whereas u+ stand for single particle pro- 



cesses 1 11 . Note that, in Eq. (H), the fields {v^(r)} are 



subject to the constraint V • v ± (r) = 0; i.e., all vortex 
lines either form closed loops or go to infinity. More im- 
portantly, it should also be noticed that the two fields 
v+ and v~ cannot be independent of each other, since 
m^(l;r) and m fl (2;r) in Eq. (Q), and hence u M (l;r) = 
[<(?) + V( f )l/2 and ^(2;F) = [«+(*) - «-(?)]/2, can 
take only integer values. As depicted with open circles 
in Fig. |^, (u+ , u~ ) at each r can take only half of the 
elements in the product set of integers ZxZ; u+ and v~ 
are topologically coupled with each other. Contrary to the 
capacitively coupled ID JJAs (llj, this topological cou- 
pling plays crucial roles in the present case, which will 
be discussed in more detail below. 

It is not difficult to understand the physics described 
by each of the Hamiltonians Hy. Unless Co = 0, 
the length-scale dependence of the anisotropy factor 
C + (q) = C(q) <C 1 is screened out at length scales larger 
than \J Ci/Cq, and thereby U^(q) is simply reduced to 
the highly anisotropic current-like interaction: U"^ (q) ~ 

Co/[A(q) + CbA(«)] and U+(q) ~ l/[A(q) + C A(uj\. 
Such an anisotropic model has been studied in Ref. JL7) , 
and is known to exhibit an anisotropic 3D transition 
which is associated with the disruption of the vortex 
loops, at K = K£ close to the 2D Berczinskii-Kosterlitz- 
Thouless (BKT) transition point |J]: K+ ~ 2/tt. In the 
case of Co = 0, the vortex lines i>+ even form 2D pancake 
vortices residing on decoupled 2D layers with U^(q) ~ 

Ci/[l + CiA(w)] andC/+(q) ~ 1/A(q)[l + CiA(w)], and 
the phase transition is precisely BKT-type. In any case, 
the system of vortex lines u+ exhibits a phase transition 
at K = K+ ~ 2/tt. On the other hand, neglecting the 
anisotropy at short-length scales, U~(r) are isotropic in 



space-time: 17,7" (q*) ^ E/ r "(q) ^ l/[A(q) + A M]- In 
consequence, it follows that the system of vortex-lines 
v~ exhibits the isotropic 3D XY-type phase transition at 
K = K~ ~ 1/2a/2. At this point, one might be tempted 
to conclude that, as K is increased, the total system Hy 
might go through two successive transitions, one at K~ 
and the other at 7Cj~ , the first of which would be ascribed 
to condensation of particle- hole pairs O]. This sce- 
nario of successive transitions, however, should be tested 
against the topological coupling discussed above between 
v+ and v~. 

For this goal, it is convenient to consider the subsystem 
{i;*} of {v+} satisfying v~ = 0. In this subsystem, v*(r) 
can take only even numbers, and hence the phase transi- 
tion could take place at K = = K+ /4 ~ 1/27T, which 
is substantially lower than K~. It means that vortices 
v* could be tightly bound even before the vortices v~ 
get bound, contradicting the assumption v~ = 0. Con- 
sequently, it follows that actual phase transition should 
take place at K c between K~ and Kj , and cannot be 
accounted for exclusively by v~ , i.e., by particle-hole 
pairs. This is distinctive different from the case of the 
capacitively coupled ID chains [O, where the vortices 
v* in analogous subsystem always form a plasma of free 
vortices regardless of K in the presumed configuration 
Co,Ci <§; 1 and the topological coupling is thus irrele- 
vant; the free vortices v* completely screen out the inter- 
action among vortices v+ . 

Nevertheless, it is evident that any corrections of the 
vortex lines v+ in the vicinity of K c is exponentially small 
in the creation energy of the smallest vortex loops (or 
nearby pancake vortex-antivortex pairs when Co = 0): 
/i+ ~ K c tt |jl8). In particular, the shift of the transition 
point K c with respect to K~ can be estimated by 



(K c -K~)/K- 



0.1, 



(11) 



where the factor 2 in the exponent is due to the topolog- 
ical coupling. 

Now I examine briefly and qualitatively the current 
drag effects in the superconducting phase, by means of 
the linear response a u ,{u>) of the current in the array I 
to the voltage applied across the array I' (see Fig. nf): 



a u ,{uj) = — lim <7#/(q,M</ 



i0 H 



(12) 



where Qw is the Fourier transform of the imaginary-time 
Green's function Q w (v,t) = (T t [I(£;t,t)I(£' ;0,0)]> 
with the time-ordered product T T and the current oper- 
ators I(£; r) = sin V x (f){£] r) (since the system is isotropic 
in x- and y-direction, only the current in the x-direction 
is considered here). Due to the symmetry between the 
two arrays, it follows that an(uj) = [a + (uj) + <t_ (u;)] /2 
and <J2\(lo) = [<J-\-(oj) — <r_(a;)] /2, where <j± are de- 
fined in a manner analogous to Eq. ( p^ ) with I^(x) = 
1(1; x) ± 1(2; x). According to the discussion above on 
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the phase transition, at K > K c , both <t+ and er_ show 
superconducting behavior: cr±{u>) = <r° (w <C 1), 

which means that the drag of supercurrents along the 
two arrays is not perfect in general. However, in the 
vicinity of the phase transition, where <j+ <C <r°_, the 
currents in two arrays can be comparable in magnitude. 
This suggests the following: The particle-hole pair is not 
so tight as in the ID case, distributing over a few lattice 
constants. Yet it is still energetically favorable enough 
to play significant (if not exclusively crucial) roles in the 
phase transition and the transport. 

Before concluding, I remark briefly on non-identical 
arrays. The difference in the intra-array capacitances 
leads to additional coupling between the vortices w+ and 
v~ with the coupling strength proportional to the dif- 
ference. The arguments on identical arrays therefore re- 



main valid qualitatively as long as 



C(l;q)+C(2;q) 



C(l;q)-C(2;q) 



The difference in Josephson coupling 

energy, on the other hand, can be effectively incorporated 
in the capacitance difference by renormalizing the param- 
eters, since all the effects considered in this work depends 
only on the relative strength of the Josephson coupling 
energy and the charging energies. 

In conclusion, quantum phase transitions in two ca- 
pacitively coupled 2D JJAs have been investigated. In 
particular, it has been found that as the coupling ca- 
pacitance increases, in appropriate parameter ranges 
{Ej/Eq,EjIEx « 1; Ej/E ,E J /E 1 « Ej/Ej < oo), 
the system exhibits an insulator-to-superconductor tran- 
sition. Contrary to the capacitively coupled ID chains, 
the transition cannot be account for exclusive by the con- 
densation of particle- hole pairs. Accordingly, the drag of 
supercurrents along the two arrays is not absolute in gen- 
eral. 
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FIG. 1. Schematic side view of the system. Each chain in 
the figure represents a 2D array. 
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FIG. 2. Topological coupling of «+ and At each 

space-time position r, (v£,v~) can take only half of the ele- 
ments in Z x Z as depicted with open circles in the figure. 
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